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Coherent electronic transport in single-molecule junctions is investigated in the Coulomb blockade 
regime. Both the transmission phase and probability are calculated for junctions with various contact 
symmetries. A dramatic suppression of the Coulomb blockade peaks is predicted for junctions where 
multiple atomic orbitals of the molecule couple to a single electrode although the charging steps are 
unaffected. 

PACS numbers: 73.23.Hk 85.35.Ds 73.63.-b, 85.65.+h, 



I. INTRODUCTION 

Coulomb blockade has been investigated in mesoscopic 
transport through metallic nanoparticlesi and semicon- 
ductor quantum dots,— and more recently in molecular 
transport 4^— Because the distance between conductance 
peaks is set by the charging energy, one might naively 
think that transport in the Coulomb blockade regime 
is essentially incoherent. Recent interferometric exper- 
iments have shown, however, that electrons transferred 
through quantum dots retain memory of their phase,^ 
a finding that has triggered intensive theoretical inves- 
tigations on the transmission phase through interacting 
systems J£r— 

Mesoscopic experiments access regimes where tunnel- 
ing transport occurs either between two discrete chan- 
nels (in lateral quantum dots) or two quasi-continua (in 
metallic nanoparticles) that are connected by a central 
quantum system with no particular symmetry. Con- 
sequently, interchannel coherence is either absent or 
washed out, and the spatial symmetry with which exter- 
nal terminals are connected to the central system matters 
little. The situation is fundamentally different in molec- 
ular transport, because molecules naturally possess sym- 
metries that strongly affect the phase coherent transmis- 
sion of electrons^ Moreover, molecular structure makes 
it possible to tunnel-couple multiple orbitals to a single 
external electrode^ 7 - 

In this article, we consider junctions composed of small 
molecules with n-fold spatial symmetry coupled to two 
external single-channel leads. The molecule-lead cou- 
pling can occur either via single or multiple spatially 
separated atomic orbitals. For single-orbital coupling we 
show that the presence of nodes in the junction's trans- 
mission spectrum depends critically upon the contact 
geometry of the junction. This is so, because the con- 
nectivity determines the relative phase accumulated by 
an electron along different, interfering paths across the 
molecule. 

The presence of such a node results in an additional 
7T shift in the junction's transmission phase d ' 11 This 
is of crucial importance for multi-orbital coupling, be- 
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FIG. 1. Transmission spectrum and average electron occu- 
pancy of a flat ethylene (C2H2) molecular junction in which 
one lead couples symmetrically to both 7r-orbitals of the 
molecule and the other couples only to one orbital. The 
HOMO resonance of the molecule yields a typical transmis- 
sion resonance, but transmission is blocked completely at the 
LUMO resonance by destructive quantum interference. De- 
spite the absence of a LUMO transmission resonance, the 
HOMO and LUMO charging steps are virtually identical. 



cause this phase shift can turn constructive interferences 
into destructive ones, thereby suppressing the conduc- 
tance. Accordingly, we find that Coulomb blockade res- 
onances can be completely suppressed in molecular junc- 
tions where several atomic orbitals of the molecule couple 
to a single electrode. 

This is the principal message of this article, and is il- 
lustrated in Fig. [1] where the transmission spectrum and 
average electronic occupancy of a flat ethylene (C2H2) 
junction is plotted against the electrochemical potential 
fi. In the flat junction, one lead couples symmetrically to 
both 7r-orbitals of the molecule while the other couples 
only to one orbital (see Fig. [5^). This junction geom- 
etry can be experimentally realized for many molecular 
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species if the molecule is deposited on a metallic sur- 
face and contacted from above by an STM tip 18 ' 19 or 
in a mechanically controlled break junction.— 7 In Fig. [TJ 
the Highest Occupied Molecular Orbital (HOMO) reso- 
nance of the molecule yields a typical transmission reso- 
nance but transmission is blocked completely at the Low- 
est Unoccupied Molecular Orbital (LUMO) resonance by 
destructive quantum interference, leading to a transmis- 
sion node. Destructive interference is not complete in 
the tails of the LUMO resonance, but the transmission is 
still suppressed by two orders of magnitude. Despite the 
absence of a LUMO transmission resonance, the LUMO 
charging step illustrated in the lower panel of Fig. Q] is es- 
sentially the same as that at the HOMO resonance. Thus 
we have the phenomenon of Coulomb blockade without a 
Coulomb blockade peak. 

This paper is organized as follows. In Sec. [H] we out- 
line the transport formalism necessary to describe single- 
molecule junctions with multi-orbital lead-molecule cou- 
pling. The linear-response transport of ethylene and 
benzene-based junctions are investigated for several junc- 
tion bonding configurations in Sec. Mil The non-linear 
transport of flat and para benzene junctions are pre- 
sented in Sec. IIVI Conclusions and discussion are pre- 
sented in Sec. [V] 



II. TRANSPORT FORMALISM 

Transport quantities are naturally expressed in terms 
of Green's functions. The retarded Green's function 
(t) = -^(^({cin^f), £^,(0)}) of the junction 
may be transformed into the energy domain giving:- 



G(E) = IE- flW - E T (#) - £c(£) 



(1) 



where H^h is the one-body part of the molecular Hamil- 
tonian H mo \, St is the tunneling self-energy matrix and 
Sc is the Coulomb self-energy matrix. H mo \ is modeled 
using a semi-empirical Hamiltonian known to accurately 
characterize the electronic spectrum of 7r-conjugated 
molecules The tunneling self-energy £t=X^ q can 
be decomposed as the sum of contributions from the var- 
ious contacts, which can be calculated exactly using an 
equations-of- motion method.^' 21 Of particular interest is 
the imaginary part of St, which causes broadening of the 
molecular resonances. Often this broadening is discussed 
in terms of the tunneling-width matrix, defined for lead 
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(2) 



The Coulomb self-energy must in general be calculated 
using an appropriate approximation. Here we utilize 
the non-perturbative method of Ref. Q to calculate Ec, 
which has been shown to accurately describe Coulomb 
blockade in single-molecule junctions. 



While the Green's function fully characterizes the 
molecular junction, the transmission phase is most read- 
ily defined in terms of the scattering matrix S a p(E). For 
purely elastic quantum transport, the dominant trans- 
port mechanism at room temperature in single-molecule 
junctions, the Green's function and scattering matrix are 
related as follows; n i 22 

S a p(E) = -5 a pl + h a {E)G{E) 1 \{E). (3) 

The off-diagonal terms of S are the transmission ampli- 
tudes [t a p(E)] nm for scattering from mode m in lead (5 
to mode n in lead a. The total transmission probability 
between leads a and /3 is given by: 



T a0 (E) = Tr{t a p(E)tl f3 (E)}, 



(4) 



where the trace is over all relevant atomic orbitals of 
the molecule. The tunneling-width amplitude matrices 
7 Q (E) are related to the tunneling-width matrices T a (E) 
as follows: 



r Q (£)= 7 t(£;) 7ct (i;). 



(5) 



In order to simplify our present discussion and analysis 
we shall henceforth: (i) Consider only two-terminal junc- 
tions, (ii) assume each lead is characterized by a single 
mode, and (iii) take the broad-band limit^i for the lead- 
molecule coupling, where r Q (£')=r Q becomes an energy- 
independent matrix of constants related to the tunneling 
rate. Together, the first two assumptions ensure the ex- 
istence of a single well-defined transmission eigenphase 
9(E)— avg [t(E)\ for the junction. 

Although the investigation of coherent quantum trans- 
port in molecular junctions is not new^ most theoretical 
studies limit their analysis to cases where T Q is diago- 
nal, corresponding physically to the special case where 
single atomic orbitals couple to single leads. Presently, 
we consider transport in a junction geometry we label 
'flat,' where a single orbital connects to one lead while 
all atomic orbitals of the molecule couple to a single chan- 
nel (e.g., s-wave) of the surface electronic states in the 
other, a configuration which is already experimentally 
accessible?! 8 ^ 



III. LINEAR-RESPONSE JUNCTION 
TRANSPORT 

In a flat junction with two leads, labeled I and II, 
the tunneling- amplitude matrices are [yi\ nm =\Ti8 na 8 ma 
and [7ii]„ m =v / rn, where a is the orbital connected to 
lead I and all orbitals couple equally to lead II. T\ and 
Tn are constants characterizing the tunneling between 
the molecule and leads I and II, respectively. Following 
Eq. ([3]), the transmission amplitude becomes: 



N 



(6) 
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FIG. 2. Numerically obtained transport results for flat, longitudinal and transverse ethylene (C2H4) junctions, (a) Schematic 
diagrams for each junction's geometry, (b) and (c) Transmission probabilities (top) and phases (bottom) for the flat (b) and 
transverse and longitudinal (c) connections, (d) Molecular 7r-orbital occupancy for all three geometries, showing that the 
destruction of the Coulomb peaks does not affect the accumulation of charge on the molecule. The transmission spectrum 
in both the longitudinal and transverse geometries is symmetric about /io, whereas it is asymmetric in the flat configuration, 
a property which should be experimentally resolvable. Using the molecular Hamiltonian and parameterization of Ref. [f|, 
the parameters for ethylene were: i=2.6eV, Un=8.9eV, Ui2=5.33eV and r=ri=rn=0.5eV in the longitudinal and transverse 
geometries while [7ii] nm =\/T/2 in the flat geometry. 



where N is the total number of molecular orbitals. From 
Eq. (jlj, we see that the transmission probability will be 
composed of diagonal terms proportional to ■ |G a _/| 2 
and cross terms proportional to J2j^k G a jG* ak . We now 
show that interferences between various transmission am- 
plitudes can cause these two terms to cancel, giving rise 
to a dramatic reduction of the transmission peaks. 

A. Ethylene junction 

For pedagogical reasons, we first consider several ethy- 
lene (C2H4) junctions, each shown schematically in 
Fig. [2^. The transmission probability and phase for the 
flat junction is shown in Fig.[5}D, while the spectra for the 
transverse and longitudinal junctions are both shown in 
Fig. In all three geometries Fi=0.5eV. In the lon- 
gitudinal and transverse cases Fn=0.5eV, while in the 
flat geometry Fn=0.25eV such that the total coupling 
f„=Tr{r n }=0.5eV. 

For chemical potentials below the center of the HOMO- 
LUMO gap (/xo), the spectra of the three geometries are 
qualitatively very similar: for every peak in one there is a 
peak in the others and similarly for the nodes, although 
the exact location of the nodes varies. For fi > /io, how- 
ever, there is a nearly 500 x reduction of the Coulomb 



blockade peak height in the flat geometry as compared 
with either the longitudinal or transverse spectra. 

The origin of this suppression can be under- 
stood by examining the transmission amplitudes in 
each geometry. In the longitudinal and trans- 
verse junctions we find t^°^ s (E)—iT [G 12(E)] and 
t^f ns (E)=iT[Gii(E)], respectively, where ri=r„=r. 
From Eq. ([5]), we see that the flat junction's ampli- 
tude tj^f (E)=iT/V2[Gu(E) +G 12 (E)] is proportional 
to the sum of the amplitudes from the other two con- 
figurations. From the longitudinal and transverse phase 
spectra, shown in the lower part of Fig. it is evident 
that the node in Gn when /i=/io causes Gn and Gi 2 to 
become it out of phase for all /i > /io. Exactly on res- 
onance with [i > /io, |Gii| = |Gi2| and consequently, via 
total destructive interference, a transmission peak in the 
longitudinal (or transverse) spectrum becomes a trans- 
mission node in the flat junction's spectrum, with a con- 
comitant lapse of the transmission phase in the latter. In 
the vicinity of a resonance, the two terms are not exactly 
equal but still interfere destructively giving rise to a pro- 
nounced reduction of the tails of the Coulomb blockade 
peaks instead of an exact cancellation. 

Transmission peaks occur at energies where a sys- 
tem has a degeneracy between two charge states. Even 
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though these Coulomb blockade peaks may be nearly de- 
stroyed by quantum interferences, as we just showed, the 
charge of the molecule still changes by one as we cross 
each resonance. The 7r-orbital molecular occupancy is 
shown for all three junction geometries as a function of 
chemical potential in Fig.[2ji and exhibits nearly the same 
spectra for all three geometries. Comparing the phase of 
transmission, shown in the bottom parts of Figs. (|2J:)15J:) 
with the charge spectra, we see that for each step in the 
occupancy there is an associated increase in the trans- 
mission phase by 7r, in agreement with the Friedel-sum 
rule^ 



B. Benzene junction 

Since the suppression of transmission peaks in fiat 
junctions are a manifestation of interference between 
multiple transport pathways, the effect should be more 
pronounced in larger molecules which possess a corre- 
spondingly larger number of transport pathways. To in- 
vestigate this hypothesis, we consider the transmission 
spectrum for a flat benzene junction in which one elec- 
trode couples to just a single 7r-orbital of the benzene 
ring, while the second couples equally to all six 7r-orbitals. 
Such a junction could be experimentally realized by de- 
positing a benzene molecule on a Pt or graphene surface 
and contacting it from above with an STM tip 18 i 19 or 
using a mechanically controlled break junction^ 

The calculated spectrum of a flat benzene junction is 
shown in Fig. [3] for several values of lead-molecule cou- 
pling. Except for the first two peaks, all the molecu- 
lar resonances coincide with transmission nodes, a fact 
which is observable from the transmission probability or 
as 7r-slips in the transmission phase, shown in the bot- 
tom portion of the same figure. Owing to the molecular 
symmetry of benzene, there are only four non-degenerate 
terms in the flat junction's transmission amplitude. Us- 
ing Eq. @: 



f Flat 



iVfifii 



in +^14^ +2 ^£12^+2^13^ 

para ortho mcta. 



(7) 



where the ortho and meta terms are doubly degener- 
ate and the dependence on energy is implicit. The first 
two peaks correspond to resonant tunneling through the 
nodeless molecular ground state, for which no destructive 
interference is possible, so that all the terms in Eq. ([7]) 
add constructively. For every other molecular resonance, 
three amplitudes are in phase and three are 7r out of 
phase such that by symmetry the total transmission am- 
plitude vanishes exactly on resonance. Elsewhere in the 
spectrum, this coordination of cancellation is visible as 
a strong suppression of the molecular resonance peak 
heights. As the lead-molecule coupling is reduced the 
individual transmission amplitudes are only appreciable 
near a resonance, where destructive interference between 




II 

w 





-71/2 
-71 



- 1 


1 








1 








1 1 




_ !? 

r 


Y 




t 












1 




-< 


1 

jj , 










, 1 





-30 -20 -10 10 20 30 
\i - il (eV) 



FIG. 3. Top panel: Sketch of the para and the flat, multi- 
orbital coupling geometries for a benzene (CeHe) molecular 
junction. Second panel: Transmission spectrum (logscale) of 
a para benzene junction. Lowest two panels: Transmission 
probability and phase for a flat benzene junction. All but the 
first two transmission peaks are completely suppressed in the 
flat junction due to destructive quantum interference. The 
destructive interference is not complete in the tails of the 
resonances, which exhibit a maximal ~3350x reduction for 
rn=0.3eV. Reductions of ~1050x and ~50x were found for 
rn=0.6eV and rn=3.0eV, respectively. The full transmission 
spectrum of the para benzene junction, shown in the top por- 
tion of this figure, is symmetric about while that of the flat 
junction is strongly asymmetric. The complete spectrum of 
each junction is included for illustrative purposes, although 
the full range of gating may not be experimentally achiev- 
able. In all cases, lead I is coupled to a single orbital with 
ri=0.6eV. In the para junction rn=0.6eV. The total coupling 
to the flat-contacted second lead Fn was varied, where each 
orbital contributed Fn/6. The model and parameterization 
used in these simulations are discussed in detail in Ref. @]. 



amplitudes is most complete, resulting in enhanced peak 
suppression with decreasing lead-molecule coupling. 



IV. NON-LINEAR JUNCTION TRANSPORT 

As we have already seen, the transmission probabil- 
ity of the junction, which determines the linear-response 
transport coefficients, is strongly affected by the con- 
tact geometry. In the following discussion, we focus 
on the non-linear junction response of a flat junction 
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FIG. 4. The non-linear current J as a function of bias volt- 
age V for para and flat benzene-Cu junctions. Because of the 
coherent suppression of the Coulomb blockade peaks in the 
fiat geometry, features in the current as a function of volt- 
age are strongly reduced as compared to the para geometry. 
A closeup of the current in the flat junction is shown in the 
top portion of the figure. In the fiat junction, the benzene 
molecule lies 2.7A above the Cu(I10) plane and 5A below the 
STM tip with lead-molecule coupling matrices rJ l I m = 1.3eV 
and rJ lm =(5„ a ^ma0.1eV, respectively, where a is an atomic 
orbital of the molecule. In the para junction both the cop- 
per surface and STM tip lie 3A from their closest respective 
atomic orbitals and each couples with I.3eV to that orbital. 
The voltage V=Vi-Vn, where Vx is the STM voltage and Vu 
is the voltage of the Cu(110) surface. 



composed of a single benzene molecule adsorbed on a 
Cu(110) surface. Although there are studies for other 
Cu surfaces^ - — we focus on the 110 plane because low- 
coverage experiments have already been performed for 
benzene using a STM setup similar to what we propose. 30 

Experimentally, it has been determined that 
A£ b <1.0eV (99kJ/mol), where AE b is the binding 
energy of benzene on Cu(llO) surfaced Using many- 
body perturbation theory the binding energy can be 
related to the lead-molecule coupling 3 -! T where, for a 
flat configuration in which the benzene molecule lies 
2.7A above the copper surface, we find Tr{r} <7.8eV, 
ca. 1.3eV per atomic orbital of the molecule. Screening 
of the intramolecular Coulomb interactions by the 
copper surface was included via the image charge 
method using a semi-empirical multipole expansion for 
the 7r-electronsi22 

The simulated non-linear responses of the para and 
flat junctions are shown in Fig. where the current 
flowing into lead a for a two-terminal device in the elastic 



cotunelling regime^ is given byM- 

/oo 
dE T a p{E) [fp(E) — f a {E)] , (8) 
-00 

where f a (E) is the Fermi distribution for lead a at chem- 
ical potential /j, a and inverse temperature /3 a . Rela- 
tive to the center of the HOMO-LUMO gap /Ucu(iio)- 
/io =_ 0.425eV, where we have used the work function of 
the Cu(110) surface^ </> Cu ( 110 )=4.48eV and ^o=(£ie + 
e EA )/2=-4.055eVi2£i 3 ii In the flat junction, the benzene 
molecule was placed 2. 7 A above the copper plane with 
the STM tip 5A above an orbital a. The tunneling-width 
matrices for the copper plane and STM tip were set to 
rJ l I m =1.3eV and rl lm =5 na S ma 0.1eV, respectively. In the 
para geometry (cf. top of Fig. [3]), the copper surface and 
STM leads were arranged opposite one another 3A away 
from their closest respective molecular orbitals with a 
tunneling-width matrix element of 1.3eV. In both junc- 
tions, the voltage was applied symmetrically. 

In the previous section, we investigated the dependence 
of the transmission function on the lead-molecule contact 
geometry. From Eq. ©, it is clear that this dependence 
directly affects the non-linear current passing through 
the junction. As shown in Fig. (Q}, the near total de- 
struction of Coulomb blockade peaks in the flat configu- 
ration gives a nearly flat I-V response, with a nearly two 
orders-of-magnitude reduction in the current step height 
as compared with the para junction. Such quantum in- 
terference effects would also affect the apparent height 37 
of adsorbed molecules in STM images. 

V. CONCLUSIONS 

In summary, we find that quantum interferences can 
effectively destroy Coulomb blockade peaks. The trans- 
mission probability of an electron tunneling through a 
flat junction, in which a single-channel lead couples to 
all n relevant atomic orbitals of a molecule with n-fold 
spatial symmetry, is determined by considering the co- 
herent superposition of transmission amplitudes from all 
possible junction connectivities. By virtue of the nodal 
structure of the many-body wavefunction, these ampli- 
tudes exactly cancel at many of the Coulomb blockade 
resonances, completely suppressing the Coulomb block- 
ade peaks but leaving the charging steps unaffected. 

Larger molecules with n-fold spatial symmetry, such as 
[18]-annulene, possesses a smaller charging energy than 
benzene or ethylene, and may therefore be preferred ex- 
perimentally for Coulomb blockade studies. We predict 
that all but the first two transmission resonances in the 
7r-electron band will be completely suppressed for a flat 
[18]-annulene junction, and that the tails of the trans- 
mission resonances will be even more strongly suppressed 
than for benzene. 

Single-molecule junctions offer a novel and interesting 
physical system in which to investigate the effects of spa- 
tial symmetry and quantum- interference on transport. 
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Since tunneling transport is central to scanning probe 
imaging, these interference effects arc likely to be impor- 
tant in interpreting images of n-fold symmetric molecules 
adsorbed on conducting surfaces. 

PJ acknowledges the support of the NSF under grant 



No. DMR-0706319. CAS thanks the Institute of Nuclear 
Theory at the University of Washington for its hospitality 
and the Department of Energy for partial support of this 
work. 



1 J. von Delft and D. Ralph, Phys. Rep. 345, 61 (2001). 

2 Y. Alhassid, Rev. Mod. Phys. 72, 895 (2000). 

3 I. L. Aleiner, P. W. Brouwer, and L. I. Glazman, Phys. 
Rep. 358, 309 (2002). 

4 S. Kubatkin, A. Danilov, M. Hjort, J. Cornil, J.-L. Bredas, 
N. Stuhr-Hansen, P. Hedegard, and T. Bjornholm, Nature 
425, 698 (2003). 

5 M. Poot, E. Osorio, K. O'Neill, J. M. Thijssen, D. Van- 
maekelbergh, C. A. van Walree, L. W. Jenneskens, and 
S. J. van der Zant, Nano Lett. 6, 1031 (2006). 

6 J. P. Bergfield and C. A. Stafford, Phys. Rev. B 79, 245125 
(2009). 

7 H. Song, Y. Kim, Y. H. Jang, H. Jeong, M. A. Reed, and 
T. Lee, Nature 462, 1039 (Dec. 2009). 

8 R. Schuster, E. Buks, M. Heiblum, D. Mahalu, V. Uman- 
sky, and H. Shtrikman, Nature 385, 417 (1997). 

9 M. Avinun-Kalish, M. Heiblum, O. Zarchin, D. Mahalu, 
and V. Umansky, Nature 436, 529 (Jul. 2005). 

10 H.-W. Lee, Phys. Rev. Lett. 82, 2358 (1999). 

11 A. Levy Yeyati and M. Biittiker, Phys. Rev. B 62, 7307 
(2000). 

12 P. G. Silvestrov and Y. Imry, Phys. Rev. Lett. 85, 2565 
(2000). 

13 G. Hackenbroich, Phys. Rep. 343, 463 (2001). 

14 A. Silva, Y. Oreg, and Y. Gefen, Phys. Rev. B 66, 195316 
(2002). 

15 C. Karrasch, T. Hecht, A. Weichselbaum, Y. Oreg, J. von 
Delft, and V. Meden, Phys. Rev. Lett. 98, 186802 (2007). 

16 D. M. Cardamone, C. A. Stafford, and S. Mazumdar, Nano 
Lett. 6, 2422 (2006). 

17 M. Kiguchi, O. Tal, S. Wohlthat, F. Pauly, M. Krieger, 
D. Djukic, J. C. Cuevas, and J. M. van Ruitenbeek, Phys. 
Rev. Lett. 101, 046801 (2008). 

18 J. Repp, G. Meyer, S. M. Stojkovic, A. Gourdon, and 
C. Joachim, Phys. Rev. Lett. 94, 026803 (2005). 

19 J. Repp, G. Meyer, S. Paavilainen, F. E. Olsson, and 



M. Persson, Science 312, 1196 (2006). 

20 C. W. M. Castleton and W. Barford, J. Chem. Phys. 117, 
3570 (2002). 

21 A.-P. Jauho, N. S. Wingreen, and Y. Meir, Phys. Rev. B 
50, 5528 (1994). 

22 D. S. Fisher and P. A. Lee, Phys. Rev. B 23, 6851 (1981). 

23 A. Nitzan and M. A. Ratner, Science 300, 1384 (2003). 

24 J. Friedel, Nuovo Cimento Suppl. 7, 287 (1958). 

25 L. J. Lauhon and W. Ho, J. Phys. Chem. A 104, 2463 
(Mar. 2000). 

26 P. S. Weiss, M. M. Kamna, T. M. Graham, and S. J. Stran- 
ick, Langmuir 14, 1284 (Mar. 1998). 

27 B. L. Rogers, J. G. Shapter, and M. J. Ford, Surface Sci- 
ence 548, 29 (Jan. 2004). 

28 A. Bilic, J. R. Reimers, N. S. Hush, R. C. Hoft, and M. J. 
Ford, J. Chem. Theory Comput. 2, 1093 (2006). 

29 H. Lesnard, N. Lorente, and M. Bocquet, J. Phys.: Con- 
dens. Matter. 20, 224012 (Jun. 2008). 

30 T. Komeda, Y. Kim, Y. Fujita, Y. Sainoo, and M. Kawai, 
J. Chem. Phys. 120, 5347 (Mar. 2004). 

31 J. P. Bergfield, J. D. Barr, and C. A. Stafford, 
arXiv:1008.0035v2. 

32 J. D. Barr, J. P. Bergfield, and C. A. Stafford, (unpub- 
lished) . 

33 M. Biittiker, Phys. Rev. Lett. 57, 1761 (1986). 

34 CRC Handbook of Chemistry and Physics, edited by D. R. 
Lide et al. (CRC Press, Boca Raton, Fla., 2005). 

35 A. Mikaya, E. Trusova, V. Zaikin, E. Karaulova, and 
L. Petrova, Bulletin of the Academy of Sciences of the 
USSR Division of Chemical Science 31, 1319 (1982). 

36 B. K. Janousek and J. I. Brauman, In: Gas Phase Ion 
Chemistry (Academic Press, New York, 1979). 

37 J. Weis, R. J. Haug, K. v. Klitzing, and K. Ploog, Phys. 
Rev. Lett. 71, 4019 (Dec 1993). 



Para 




2x 



lat 



-30 -20 -10 10 
H - n (eV) 



20 



'i 



30 



